Prom target to projectile and back again: selfduality of high energy evolution. 
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We prove that the complete kernel for the high energy evolution in QCD must be selfdual. The 
relevant duality transformation is formulated in precise mathematical terms and is shown to trans- 
form the charge density into the functional derivative with respect to the single-gluon scattering 
matrix. This transformation interchanges the high and the low density regimes. We demostrate 
that the original JIMWLK kernel, valid at large density is indeed dual to the low denisity limit of 
^ ■ the complete kernel derived recently in hep-ph/0501198 
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Recent months have seen renewed attempts to understand Pomeron loops contribution to the high energy evolution 
of hadronic cross section in QCD 0. Over the years much work has been done to understand the evolution in the 
reggeon language 0- In recent years the study of the high energy scattering has centered around the so called JIMWLK 
evolution equation [Id LI and its mean field limit 0. These describe the approach of the scattering amplitude to 
saturation due to multiple scattering corrections on dense hadronic targets, or in the diagrammatic language, the fan 
diagrams. The JIMWLK equation however only partially takes into account the processes whereby the gluons emitted 
in the projectile wave function at an early stage of the evolution, are "bleached" by subsequently emitted gluons, or 
. the so-called Pomeron loops. 

The work 0, lead to realization of the absence of the Pomeron loops in the JIMWLK evolution equation and 
the resulting violation of i-channel unitarity. Very recently several papers have derived corrections to the JIMWLK 
evolution equation away from the high density limit 0, E| • The most complete of the results so far is given in [ljj 
\ where the kernel for the evolution equation in the low density limit is derived. This kernel has a very specific form 
■ which suggests an intriguing duality between the high density and low density limit. 

In the present paper we prove that such a duality of the evolution kernel is in fact required by the i-channel 
unitarity and the Lorentz invariance of the theory. We consider the evolution from the projectile side and the 
' target side separately. The requirement that the evolution of the projectile and the target wave functions has the 
same functional form coupled with the requirement of Lorentz invariance of the scattering matrix, leads to the 
f-} • condition that the kernel of the evolution x\Pt j~] must be invariant under the transformation 
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where p is the charge density in the target wave function, and a is the phase of the scattering matrix of a single gluon 



x. 

on the same target (see below). The discussion of this paper assumes the eikonal form of the S- matrix, and is thus 
valid in the leading logarithmic approximation. 

We start by considering a general expression for the S'-matrix of a projectile with the wave function \P) on a target 
with the wave function |T)|l2(, where the total rapidity of the process is Y. The projectile is assumed to be moving 
to the left with total rapidity Y — Yq (and thus has sizeable color charge density p~), while the target is moving to the 
right with total rapidity Yq (and has large p + ). We assume that the projectile and the target contain only partons 
with large k~ and k + momenta respectively: k~ > A~ and k + > A + . The eikonal expression for the S'-matrix reads 



S Y = J Dp+ a (x,x-) W$ \p + {x-,x)} ££_y [a], (2) 
where S p is the S'-matrix averaged over the projectile wave function|l4| 

(V denotes path ordering with respect to x~), and W T [a] is the weight function representing the target, which is 
related to the target wave function in the following way: for an arbitrary operator O[y0 + ] 

(T\6[p+(x)}\T) = [ Dp+ a W T [p+(x-,x)} 0[p+(x,x-)}. (4) 



The field a(x) is the A + component of the vector potential in the light cone gauge A — 0. This is the natural gauge 
from the point of view of partonic interpretation of the projectile wave function. In the same gauge the target charge 
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density p + is related to a through the solution of the classical equations of motion 
a a (x,x-)T a = ±(x-y) {U\y,x-) p +a (y,x-)T a U(y,x-)} , U{x,x~) = V exp{t j dy~T a a a (x,y-)} . 



(5) 

where T b a c = if abc is the generator of the SU(N) group in the adjoint representation. The unitary matrix U(x,x — > 
oo ) has the meaning of the scattering matrix of a single gluon from the projectile wave function on the target \T). 
In the formulae above we use hats to denote quantum operators. Note that the quantum operators 

p- a (x)= dk-al b (x,k-)T£ c a?(x,k-), p +a (x) = / dk+af{x, k + )T£ c a c t (x, k+) (6) 
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do not depend on longitudinal coordinates, but only on transverse coordinates x (here a* and a* stand for single gluon 
creation operators of left and right moving gluons respectively). The "classical" variables a and p + on the other hand 
do depend on the longitudinal coordinate x~ . This dependence, as discussed in detail in arises due to the need 
to take correctly into account the proper ordering of noncommuting quantum operators. Thus the ordering of the 
quantum operators p + in the expansion of O in the lhs of eq. 10} translates into the same ordering with respect to the 
longitudinal coordinate x~ of p + (x~) in the expansion of 0[p + (x~)} in the rhs of eq.@. 

As was shown in the functional W T [a] cannot in general be interpreted as probability density, as it contains a 
complex factor. This factor - the Wess-Zumino term, ensures correct commutators between the quantum operators 
p a . In the present paper we do not require an explicit form of this term, but the following property which is implicit in 
eq. (J3J is crucial to our discussion. The "correlators" of the charge density (p ai {x\, xV)...p a " (x n , x~ )) do not depend 
on the values of the longitudinal coordinates x~ , but only on their ordering jlOj. 

Note that one can define an analog of W T for the wave function of the projectile via 

(P\6[p-{x)]\P) = J Dp- a W p [p-]0[p-{x,x-)]. (7) 

With this definition it is straightforward to see that E p and W p are related through a functional Fourier transform. 
To represent E as a functional integral with weight W p we have to order the factors of the charge density p~ in the 
expansion of eq.©, and then endow the charge density p~ (x) with an additional coordinate t to turn it into a classical 
variable. This task is made easy by the fact that the ordering of p in eq.© follows automatically the ordering of the 
coordinate x~ in the path ordered exponential. Since the correlators of p(x, ti) with the weight W p depend only on 
the ordering of the coordinates ti and not their values, we can simply set t = x~ . Once we have turned the quantum 
operators p into the classical variables p(x~), the path ordering plays no role anymore, and we thus have 

E p (a) = j Dp a W p [p] e «/«to~/««V"(*.0« a (*.«~). ( 8 ) 

We now turn to the discussion of the evolution. The evolution to higher energy can be achieved by boosting either 
the projectile or the target. The resulting S-matrix should be the same. This is required by the Lorentz invariance of 
the iS-matrix. Consider first boosting the projectile by a small rapidity SY. This transformation leads to the change 
of the projectile S'-matrix E of the form 

A E P =x t [a ,^ ]E P [a] (9) 

The explicit form of the kernel \ 1S known in the large density limit 3,0 and in the small density limit [To| . 
Substituting eq.© into eq.JSJ we have 
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Dp+ a (x,x-) Wl\p + {x-,x)] {x f K^] S£-Y M 



Dp+ a (x,x-) \ x [ a ,^-}W p [ P + (x-,x)}\ E£_ y >]. (10) 



Where the second equality follows by integration by parts. We now impose the requirement that the S'-matrix does 
not depend on Y [13|, Since E in eq.© depends on the difference of rapidities, eq.© implies 

A-E = - X t K i. ]sH . 
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Requiring that dS/dYo — we find that W should satisfy 



(12) 



Thus we have determined the evolution of the target ea. l|12ll by boosting the projectile and requiring Lorentz 
invariance of the S'-matrix. On the other hand the extra energy due to boost can be deposited in the target rather 
than in the projectile. How does W T change under boost of the target wave function? To answer this question we 
consider the relation between E and W together with the evolution of E. Referring to eqs.JSJ) and (0 it is obvious 
that multiplication of E p by a is equivalent to acting on W p by the operator —iS/5p, and acting on E p by S /5a is 
equivalent to multiplying W p by ip. Additionally, the action of ip and —i5/5p on W p must be in the reverse order 
to the action of 5/ 5a and a on E p . This means that the evolution of the functional W p is given by 
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W p = xl-ij-M W p [p\. 
op 



(13) 



Although ca. l|13(l refers to the weight functional representing the projectile wave function, the i-channel unitarity 
requires that the functional form of the evolution must be the same for W T . Comparing ea. l|12|) and ea. (|13|l we find 
that the high energy evolution kernel must, as advertised, satisfy the self duality relation 



(14) 



This is the main result of the present paper. 

It is interesting that although we do not know the complete kernel Xj we can check that the known limits of it are 
consistent with the duality eq.QJ. In particular, the duality transformation interchanges the large and small density 
limits, and the kernel \ is known in both these limits. In the large density limit the pertinent expression is the 
JIMWLK kernel. It can be conveniently written as ^3 
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The low density limit of the kernel has recently been calculated in 10] . It reads 
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Obviously the kernel ea. (|15|l transforms into the kernel ea. l|16|l by the duality transformation eq.Q. 

One can go one step further and consider also preasymptotic terms. In [lj| we have derived a set of preasymtotic 
corrections to ea. (115(1 (not a complete set) which take into account higher density terms in the projectile wave function: 
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The field bf satisfies the "classical" equation of motion 

{d l + 9f ahc K{z, [p])} d+b^z, [p]) - gp a (z, x-) ; 

e^dib^z, [p]) - djbUz, [p]) + gf ahc bt(z, \p])b](z, [p])] = . (18) 

A correction to the low density limit was derived in [Toj]. resumming the same type of corrections but in the target 
wave function. The relevant expression is 

1 
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Again, we see that the two expressions eas. (|17p > and (|19|l are related by the duality transformation eq.Q). 

In this Letter we have established the exact selfduality property that has to be satisfied by the complete kernel 
of the high energy evolution in the leading logarithmic approximation. The existing expressions of the limits of \ 
are transformed into each other by the duality transformation, as they should. We hope that this property will be 
helpful in future work aimed at finding the complete kernel. The selfduality must be preserved by any expression 
which purports to contain all Pomeron loops. In fact, one can suggest an expression based on the known large and 
small density asymptotics: 
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This provides an interpolation between the large and small field limits, and preserves selfduality. Of course, this 
interpolation is not unique, and the exact kernel may turn out to be different. The expression eq. l|2UI) has however 
a simple property, that it becomes the original JIMWLK expression rotated into a different basis. Physically the 
meaning of the interpolation ea. l|2(J|) can be understood by recalling that iS/5a has the meaning of p~ [l^. Thus 
eq. (|20|l suggests that the physics does not depend on the current densities p + and p~ separately, but only on the total 
color charge density in the lab frame p° = p~ + p + . 

The selfduality of the kernel is somewhat similar (although different in detail) to the duality symmetry of a harmonic 
oscillator Hamiltonian p —> x, x — > —p. One thus hopes that it may eventually be of help in solving the complete 
evolution equation, once it is derived. 

Finally we would like to mention that Lipatov's effective action approach 15] to the high energy evolution in QCD 
is claimed to respect Lorentz invariance and both the direct and cross channel unitarity. It would be interesting to 
relate Lipatov's effective action as well as reggeon field theory directly to our present approach. 
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